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The multipole moments and the tidal Love numbers of neutron stars and quark stars satisfy 
certain relations which are almost insensitive to the star’s internal structure. A natural question is 
whether the same relations hold for different compact objects and how they possibly approach the 
black-hole limit. Here we consider “gravastars”, which are hypothetical compact objects sustained 
by their internal vacuum energy. Such solutions have been proposed as exotic alternatives to the 
black-hole paradigm because they can be as compact as black holes and exist in any mass range. By 
constructing slowly-rotating, thin-shell gravastars to quadratic order in the spin, we compute the 
moment of inertia /, the mass quadrupole moment Q, and the tidal Love number A in exact form. 

The I-X-Q relations of a gravastar are dramatically different from those of an ordinary compact 
star, but the black-hole limit is continuous, i.e. these quantities approach their Kerr counterparts 
when the compactness is maximum. Therefore, such relations can be used to discern a gravastar 
from an ordinary compact star, but not to break the degeneracy with the black-hole case. Based 
on these results, we conjecture that the full multipolar structure and the tidal deformability of 
a spinning, ultracompact gravastar are identical to those of a Kerr black hole. The approach to 
the black-hole limit is nonpolynomial, thus differing from the critical behavior recently found for 
strongly anisotropic neutron stars. 

PACS numbers: 04.20.-q, 04.25.-g, 04.70.Bw, 04.30.-w. 


I. INTRODUCTION 

Astrophysical black holes (BHs) are the simplest 
macroscopic objects in the Universe, being characterized 
only by their mass M and angular momentum J. This 
property is formally proved by the no-hair and unique¬ 
ness theorems Erlll which - roughly speaking - state 
that any regular BH cannot possess further asymptotic 
charges other than M and J, and that any stationary, 
vacuum solution of Einstein’s equations is described by 
the Kerr metric. This “two-hair” property implies that 
all multipole moments of a Kerr BH can be written in 
terms of M and J in an elegant and compact form Q. 
The simplicity of BHs survives also in the presence of 
external tidal fields. The tidal Love numbers (which 
measure the tidal deformability of a self-gravitating ob¬ 
ject [3, 0) of nonspinning [fj-llH and spinning [l2Hl3 | 
BHs are precisely zercf] 

It has been recently discovered that relativistic com¬ 
pact stars satisfy certain nearly-universal relations which 
are reminiscent of the no-hair properties of BHs. The 
multipole moments and the tidal Love numbers of a 
neutron star (NS) depend on each other through rela¬ 
tions which are almost insensitive to the star’s internal 
structure, i.e. they depend only mildly (at the percent 
level) on the equation of state of matter in the NS inte¬ 
rior [k| [l6[ (cf. also Refs. [I3~Qjl for some earlier related 
work). 


* paolo.pani@romal.infn.it 

1 At least in the axisymmetric case to second order in the spin El 
and generically to first order in the spin E3- 


In their simplest form these relations involve the mo¬ 
ment of inertia /, the mass quadrupole moment Q , and 
the electric, quadrupolar tidal Love number A. For this 
reason, they were dubbed “I-Love-Q” relations mm. 
Although the I-Love-Q relations were originally discov¬ 
ered for slowly-rotating, barotropic, isotropic, unmagne¬ 
tized and isolated stars [15, 16], they have been extended 
to include rapid rotation mu, nonbarotropic 0] and 
| fluids, strong magnetic fields p3 |. 


anisotropic [25l 
dynamical configurations m and also deviations from 
General Relativity 0, Ebl . 2S - l31 |. The outcome of these 
studies is that the approximate universality is remarkably 
robust in realistic configurations. 


It is therefore important to explain the origin of these 
properties. Some possible explanations were proposed 
in [32l , ]33j in terms of an emergent approximate self- 
similarity in isodensity contours inside the star and in 
terms of the approach to the incompressible limit of re¬ 
alistic equations of state, respectively. Some analytical 
models to support these arguments were developed in 
Refs. 0, .35] and Ref. 0, respectively. 

Two natural questions arising in this context are 
whether the same nearly-universal relations exist for 
compact objects other than NSs and quark stars, and 
whether such approximate universality is related to the 
approach to the BH limit for very compact objects. In 
this paper we wish to consider these two problems. 


A major obstacle in understanding the approach to the 
BH limit consists in the fact that ordinary compact stars 
have a maximum mass and a maximum compactness and, 
therefore, cannot be as compact as a BH. In order to cir¬ 
cumvent this problem, Refs. 0, 0 recently studied the 
I-Love-Q relations for anisotropic, incompressible stars, 
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which can almost reach the BH limit for large values of 
the anisotropic parameter. They found that the approach 
to the BH limit is continuous and reminiscent of a phase 
transition, with the multipole moments displaying a cer¬ 
tain critical behavior. 

In this paper, we study the I-Love-Q relations for a dif¬ 
ferent model of compact object, namely for “gravitational 
condensate stars”, or gravastars [13] • In these models, 
the spacetime is assumed to undergo a quantum phase 
transition in the vicinity of the would-be event horizon. 
The horizon is effectively replaced by a transition layer 
and the BH interior is replaced by a patch of de Sitter 
space [38]. The effective negative pressure of the de Sit¬ 
ter interior contributes to sustain the self-gravity of the 
object for any compactness. In the static case these mod¬ 
els are thermodynamically [13] and dynamically (39l44l[ 
stable to linear ordeiQ for reasonable equations of state 
of the transition layer. 

In the context of the approximate universal relations 
mentioned above, gravastars have various interesting fea¬ 
tures. First, these solutions can exist for any compact¬ 
ness, smoothly connecting the Newtonian regime to the 
BH limit. As we will show, the model is sufficiently sim¬ 
ple to admit an exact solution, so that the I-Love-Q rela¬ 
tions for a gravastar can be computed analytically. This 
also sheds some light on the approach to the BH limit, 
which is not polynomial in this model but involves log¬ 
arithmic terms. Secondly, in general there is no reason 
to expect that the BH limit of a compact material body 
be continuous and one could hope to use the I-Love-Q 
relations to break the degeneracy between BHs and “BH 
mimickers”. We show that this is not the case for gravas¬ 
tars, because their multipolar structure and their tidal 
deformability in the BH limit are precisely the same as 
those of a Kerr BH. 

Through this work we use geometrized units where 
Newton’s constant and the speed of light are set to unity; 
we use c := M/R to denote the compactness of the ob¬ 
ject. 


II. SLOWLY-ROTATING THIN-SHELL 
GRAVASTARS 

We focus on the simplest gravastar model HU, which 
consists of an exterior, vacuum solution of Einstein’s 
equation glued to an interior solution supported by a 
positive cosmological constant (see also Refs. [33, 113113 
for different models). The interface defines the object’s 
radius R; a material thin shell located at r = R is re¬ 
quired because the matching of the two geometries is not 
smooth 1353. 


2 The arguments presented in Refs. 42. [fV suggest that ultracom- 

pact objects without an event horizon might be nonlinearly un¬ 
stable. Furthermore, if highly spinning these objects are linearly 
unstable against the ergoregion instability 144 . |45|| (cf. Ref. I4q| 
for a review). 


Following Hartle and Thorne [49j, [50] , we consider the 
following ansatz for a slowly-rotating object to second 
order in the spin, 

Y I 2(™o-t-T7»2-P2) 

ds 2 = -/ [1 + 2 (jo + j 2 P 2 )] dt 2 +-- dr 2 

+r 2 [1 + 2(i> 2 - j 2 )P 2 ] [dd 2 + sin 2 d(d<p - utdt) 2 ] , (1) 

where P 2 := P 2 (cost?) = (3cos 2 t? — l)/2 is a Legendre 
polynomial. The radial functions / and h are of zeroth 
order in rotation, oj is of first order, and jo, j 2 , mo, m 2 , 
V 2 are of second order. 

Using the ansatz above, it is straightforward to solve 
Einstein’s equations perturbatively for x := J/M 2 -C 
1. To zeroth order in the spin, Birkhoff’s theorem 
guarantees that the exterior solution is identical to the 
Schwarzschild metric, whereas the interior is described 
by a patch of de Sitter spacetime. Thus, to O(x 0 ), the 
solution reads furl [ffil 


f(r) = h(r) 


1 — 2cr 2 /R 2 r < R 
1 - 2 M/r r > R ’ ^ 


where M is the gravastar mass measured by an observer 
at infinity, c := M/R is the compactness, and the ef¬ 
fective cosmological constant of the de Sitter region is 
A := 8np\ = 6 M/R 3 . The junction conditions at the 
surface (r = R) have already been partially chosen by 
requiring the induced metric <jij (where Li = t,d,p) to 
be continuous across the shell (cf. Ref. [4d| for details). 
Israel’s junction conditions then relate the disconti¬ 
nuities of the extrinsic curvature on the shell with the 
stress-energy tensor of the thin layer. From these condi¬ 
tions, the surface energy E and surface tension 0 of the 
shell read [HJ 


[[VK]] = -4ttPE , 


f'Vh/f 


8tt(E - 20), (3) 


where the symbol ][A(r)]] := lim e _s. 0 [A(P + e) — A(R — 
e)] denotes the discontinuity of a generic function A(r) 
across the shell. In this simple model, the coefficient 
h is continuous across the shell, and therefore E = 0, 
whereas the surface tension is nonzero, 0 = — - r ,-’/ 

87rR\/l—2c 

It is easy to see that the thin-shcll matter satisfies the 
weak energy condition but violates the dominant energy 
condition. Furthermore, the surface tension diverges in 
the BH limit, as expected by the fact that in this case 
an event horizon is located at r = 2 M and the future 
domain of dependence of the interior has no intersection 
with the exterior. 

Note that the function / can in general be discontin¬ 
uous if the time coordinate of the external patch is dif¬ 
ferent from the one in the internal patch. This freedom 
corresponds to different choices of the equation of state 
of the thin shell. However, for all choices different from 
the one we made, the corresponding solution has a max¬ 
imum compactness smaller than that of a Schwarzschild 
BH (cf. e.g. Ref. [52|). Since in this paper we are inter¬ 
ested in the approach to the BH limit, we will impose that 
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the time coordinate is the same in the two patches, which 
singles out the specific equation of state presented above, 
allowing for models with c < 1/2. A detailed analysis of 
different equations of state will appear elsewhere [531 ] . 

At variance with ordinary stars, gravastars are purely 
gravitational solutions which are not supported by the 
pressure of internal fluids. Rotation does not induce any 
fluid motion and it is thus particularly easy to obtain 
the first- and second-order spin corrections to the static 
solution above. To first-order in the spin, the gravito- 
magnetic term satisfies the following ordinary differential 
equation 


w" + -w' = 0 , (4) 

r 

both in the interior and in the exterior, where a prime 
denotes derivative with respect to r. The general solu¬ 
tion to this equation is w(r) = cq + 2ci/r 3 , where c, are 
integration constants. One can fix one of these constants 
both in the interior and in the exterior by requiring reg¬ 
ularity at the center and asymptotic flatness at infinity, 
respectively. We obtain 


w(r) 


n r < R 

yir r > R 


(5) 

I 


where f1 := cq and J := ci are now the object’s angular 
velocity and angular momentum, respectively, and are 
still free constants. However, continuity of the induced 
metric across the shell implies ui to be also continuous, 
and this requirement translates into a relation between 
f l and J, namely f 1 = 2 J/R 3 . Therefore, to 0(x) the 
moment of inertia of a gravastar is simply 


J _ R 3 

n ~ 2 ’ 


( 6 ) 


and reduces to the Kerr value, I = 4M 3 + 0(y 2 ), when 
c= 1/2. 

Let us now compute the 0(x 2 ) corrections. In this 
case it is more convenient to discuss the interior and the 
exterior solutions separately. The exterior geometry is 
the generic (stationary, axisymmetric) vacuum solution 
of Einstein’s equation to second order in the spin, first 
derived in Refs. [H,[5C|, namely 


T r 4 /? 4 

mo = X — —— 

m 2 = ^ ^ Sq ^ cR(cR — r) (3 r 2 — 6crR — 2 c 2 R 2 ) + 3r 2 (r — 2 cR) 2 tanlW 1 
3 (r - 5cR)(r - 2cR) 


2 c 2 rR 2 
2 c 5 R 5 


cR 


r — cR 


c a R 4 — r 3 Sm. 


Jo = X — 


r 3 (r — 2 cR) 
2 


J 2 = 


X 


2c 2 r A R 2 
2 


2 c b R b [r + cR) + 


3 Sq (cR(cR — r) (6 crR + 2 c 2 R 2 — 3r 2 ) — 3r 2 (r — 2 cR) 2 tanh 1 ) 


V2 = 


2L 

cR 


Sq 3 (r — cR) tanh 


-l 


cR 


r — cR 


r — 2cR 

cR (3r 2 — 6 crR + c 2 R 2 ) \ c 5 R 5 

r(r — 2 cR) 1 r 4 


(7) 

( 8 ) 

(9) 

,( 10 ) 

( 11 ) 


where Sm and Sq are dimensionless integration constants 
related to the spin-induced shifts of the mass and of the 
quadrupole moment of the object, respectively [49|, [50}. 
In the interior, r < R, the spacetime is a solution of 


Einstein’s equation with a positive cosmological constant, 
A = 6c/ R 2 . The solution that is regular at the center 
reads 





























4 


mo = 0 

m 2 = - (^V^rR (lOcr 2 - 3 R 2 ) +3 (R 2 - 2cr 2 ) 2 tanh' 1 

jo = c 4 i? 4 x 2 a 0 , 

j 2 =- -J-—^ - | V2crR (lOcr 2 — 3i? 2 ) + 3 (i? 2 — 2cr 2 ) 2 tanh” 1 

96/2r 3 (2cr 2 - i? 2 ) ^ 1 7 V ' 

/2cr 


<22 , 


a/ 2cr 


i? 


&2 , 


t <2 =- 7 =r——^- | /2cr (4cr 2 — 3i? 2 ) + (3i? 3 — 6 cr 2 R) tanh 1 

24/2(2 cr 3 -rR 2 )\ V ' V ' 


i? 


a 2 , 


( 12 ) 

(13) 

(14) 

(15) 

(16) 


where ao and a 2 are integration constants. 

Continuity of the induced metric gij across the shell re- 

I 


lates the integration constants of the two solutions above, 
namely: 


c 4 R — 5m 

a ° = c 4 (l - 2 c)R 5 ’ 
96c 3 / 2 


-1 


1 — c 


where 


a 2 = ~ ) (^ c (l^c 2 — 3c — 3) + 3 (l — 6c 2 + 4c 3 ) tanh 

5 

5q = ^ ^ ^2 y/c (16c 2 — 9 — 6c) — 9/2 (4c 2 — l) tanh” 1 /2c ^ , 


A(c) = 2/H ( c (9 - 12c + 9c 2 + 8c 3 ) - 3 (3 - 7c + 6c 2 ) tanh' 


1 — c 


—3/2 tanh 1 /2c f c (3 — 6c — 5c 2 + 6c 3 ) — 3 (l — 3c + 4c 3 ) tanh 1 


1 — c 


(17) 

(18) 
(19) 


( 20 ) 


As expected, matching the interior geometry with the 
exterior geometry yields a solution which is continuous 
across the shell, except for the g rr component whose dis¬ 
continuity is related to the properties of the shell through 
the junction conditions. The induced metric is continu¬ 
ous but its derivatives are generically discontinuous, the 
jumps across the shell being related to the physical prop¬ 
erties of the latter. Once a model for the thin shell is as¬ 
sumed, its properties and the parameter 5m are uniquely 
determined in terms of the discontinuities of the extrin¬ 
sic curvature (3Sj, [Hl|. A more detailed discussion of the 
properties of the thin shell for various equations of state 
is in preparation IaIFI. 

The spin-induced mass quadrupole moment Q of the 
solution can be extracted from the large-distance behav¬ 
ior of j 2l namely j 2 —> Q/r 3 as r —> 00 [49|, [50j. From 


3 Very recently, Ref. [H3 l investigated in detail the junction 

conditions for a slowly-rotating, thin-shell gravastar (see also 
Refs. I 54 I . [55l for some related studies on rotating thin shells 

in general relativity.) 


Eq. (flOl) . we obtain 


Q(c) 



( 21 ) 


where 5q(c) is given in Eqs. (fTiJI) (l20l) . 

It is straightforward to evaluate the quadrupole mo¬ 
ment (1211) in the Newtonian (c —> 0) and the BH (c —> 
1/2) limits, i.e. 


fS-Wc+Wc 2 c ^° 

l -§[log(l-2c)] 1 1/2 


( 22 ) 


Note that Q is a monotonic function of c: it is negative 
(i.e. the gravastar is prolate) in the Newtonian limit, 
vanishes at c k 0.358, and becomes positive (i.e. the 
gravastar is oblate) for larger compactness. Finally, 5q —► 
0 and Q —> y 2 M 3 in the BH limit, precisely as in the 
Kerr case. This property seems to be valid also for more 
generic equations of state of the thin shell, as shown in a 
recent work (5 J . 
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III. TIDAL PERTURBATIONS OF A 
THIN-SHELL GRAVASTAR 


The tidal deformability of a self-gravitating object is 
encoded in the tidal Love numbers, which measure the 
deformation of the multipolc moments of the object in¬ 
duced by an external tidal field (see e.g. [l, jjj). For 
relativistic, nonspinning compact objects, the tidal Love 
numbers can be computed through the method devel¬ 
oped in Ref. [5(| (see also Refs. [H [lOj])- Here we fo¬ 
cus the electric tidal Love numbers of a static gravastar; 
the extensions to the magnetic SHHH and slowly- 
spinning |T3 cases are straightforward. 

We consider static, polar perturbations of the nonro¬ 
tating background defined by the metric © with x = 0 
(i.e. where all functions other than / = h vanish). 
Following Ref. j5(|, we expand the metric as g a b = 

g ( ab + s 9ab, where (a, b = t,r, d, ip), g^ is the background 
geometry anq3 

Sg ab = diag (^fH^ e \r), h~ l H [ 2 n {r), 

r 2 K w (r),r 2 sin 2 0KW{r))YtQ(0,<p). (23) 


In the above equation Y^ m are the spherical harmonics 
and a sum over £ is implicit. Due to the spherical sym¬ 
metry of the background, the azimuthal number m is 
degenerate and can therefore be set to zero without loss 
of generality. Furthermore, multipoles with different val¬ 
ues of £ decouple. By plugging the perturbation above 
into the linearized field equations, we obtain H 2 = H q 
am0 


K^(r) = 


36 M 2 


3 {£{£ + 1) - 2) r (§M - 3r + r 3 A) 

+18 (£(£ + 1) - 1) Mr - 9 (£(£ + 1) - 2) r 2 5 
+3r 3 (l6M + (£{£ + 1) - 4) r) A - 2r 6 A 2 ) H ( 


W 


-2 r (r 3 A - 3m) (&M - 3r + r 3 A) H { 


) 


rW 


, (24) 


whereas H, 


W 


satisfies a second-order differential equation 


d 2 H^ 

dr 2 


- VH& = 0 , 


(25) 


4 Similarly to the background solution, here we assume that the 
time coordinates in the interior and in the exterior of the per¬ 
turbed solution are the same. This effectively corresponds to 
fixing the speed of sound on the shell. Other choices will be 
discussed elsewhere [53ll . 

5 In order to treat the interior and the exterior solutions at the 

same time, we have derived the equations for the background 
quantities f = h = 1 — 2 M/r — Ar 2 /3, where we introduced the 
definitions A = A 0(R — r ) and M = M©(r — R ), ©(#) being 
the Heaviside step function. Therefore, A = 0 when r > R and 
M = 0 when r < R. In both cases the equations presented in 
the main text simplify considerably. 


where 


R 4 


£{£ + 1)/ + (9 - r 2 Aj 


, (26) 


and the tortoise coordinate is defined such that dr/dr* := 
g{ r ) = -^2 ^6 M + Ar 3 — 3r'j. Note that the potential V 

is discontinuous across the radius, whereas g(r) is con¬ 
tinuous by virtue of the fact that A = 6 M/R 3 . 

Because of such discontinuity, appropriate junction 
conditions have to be derived for the function Hq . Since 
Eq. (1251) is homogeneous and linear, we can assume that 
Hq is continuous across the shell without loss of gener- 

ality. To obtain the discontinuity of ILq at the radius, 
we integrate Eq. (l?5l) from R — e and R + e and then 
take the e —> 0 limit. This procedure yields the junction 
condition 


H, 




R 

9(R) 




(27) 


where we recall that a prime denotes derivative with re¬ 
spect to r. The discontinuity of the potential does not 
depend on £, [[V]] = 108 M(M — R)/R 4 , and therefore 



36 


c(l — C )H«\R) 
1-2 c R 


(28) 


for any value of £. 

Nonradial perturbations of a static, thin-shell gravas¬ 
tar have been studied in detail in Refs. (4(1 4lJ, 4J, j59|. 
The perturbation equations in the interior can be solved 
analytically in terms of hypergeometric functions, even 
in the time-dependent case [U, [43!, H§|- Hr the static 
limit, an analytical solution is also available in the ex¬ 
terior, so that the problem of static tidal perturbations 
of a gravastar can be solved in exact form. Indeed, it is 
easy to show that the regular solution of Eq. (l25l) reads 


m = 


Oil 


( rWA\ f 


=4 4 j 
2 ' 2' 1 
2/ 


+ 


2 ’ 3 


a 2 P 2 (^-l)+/? 2 Q 2 (^-l) 


r < R 
r > R 


(29) 

where P 2 and Q 2 are associated Legendre functions, 2 Fi 
is the hypergeometric function, and and /3j are inte¬ 
gration constants. One integration constant associated 
to the inner solution has been already fixed by requiring 
regularity at the center [52 ■ Other two constants (say 

= 0 and 


W 


a 2 and /3 2 ) can be fixed by imposing H ( 

Eq. (|28l) . thus leaving a solution which is defined mod¬ 
ulo an overall normalization factor, as expected from the 
linearity and homogeneity of Eq. (1251) . 

The solution above is valid for generic values of £. For 
simplicity, let us focus on the dominant quadrupolar case, 
fixing 1 = 2. At large distances, the solution behaves as 


H, 


( 2 ) 


3ci w- 6n n 


8 M 3 


c 2 + 0((A//r) 4 ), (30) 
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where Ci and C 2 are constants related to ai and /? 2 - After 
the matching discussed above, these constants are both 
proportional to the overall normalization factor. Using 
the same definition adopted in Ref. [56j, we define the 


electric, quadrupolar tidal Love number as 

8 c 2 

5 ci 


A=°^M 5 . 


(31) 


Note that A depends on the ratio C 2 /C 1 , and therefore the 
normalization factor cancels out, as expected. Finally, by 


imposing 


H, 


W 


= 0 and Eq. (l28l) . we can explicitly 


compute the value of C 2 /C 1 , yielding 


_ 56(1-2c) 2 
A - r(c) 


lOy^ (72c 3 


108c 2 + 34c - 3) + 3v / 2(l - 2c) 2 (36c 2 - 40c + 5) tanh” 1 y/2c 


(32) 


where 

r(c) = 210cVc(720c 6 + 440c 5 - 5084c 4 + 6590c 3 - 3210c 2 + 645c - 45) + 315 ( X ~ 2c ) 

x |lOVc(6c - l)(12c 2 - 16c + 3) log(l - 2c) + v^tanh” 1 y/2 c [2c(72c 5 + 72c 4 - 474c 3 + 488c 2 - 165c + 15) 
+3(1 - 2c) 2 (36c 2 - 40c + 5) log(l - 2c)] } . (33) 


The result above is exact and valid for any compactness. 
In the Newtonian and BH limits the tidal Love number 
reduces to 

A f c —> 0 

^^\_i|8( c _i) 2 c —>■ 1/2 ' (34) 

In the Newtonian limit the scaling of the tidal Love num¬ 
ber with the compactness is different than in the stellar 
case, where A/M 5 ~ c - 5 0E1. Furthermore, at vari¬ 
ance with the NS case, the electric tidal Love number 
is negative for any compactness (i.e. an oblate exter¬ 
nal quadrupolar tidal field makes a spherically-symmetric 
gravastar prolate and vice versa), and vanishes in the BH 
limit, i.e. A —> 0 as c —>• 1/2. 

These two peculiar properties are captured by a simple 
Newtonian model. For a self-gravitating barotropic fluid 
in the Newtonian regime, the quadrupolar electric Love 
number is governed by the equation J56| 


H, 


( 2 )' 


-Hf 

r 


47 rp 


p'(p) 


h^ 2) = 0 . 


(35) 


Although in dynamical situations the interior de Sitter 
spacetime is not equivalent to a barotropic fluid, one 
can attempt at modeling the gravitational perturbations 
through the equation of state P(p) = —p = — A/(87r). In 
such case Eq. ( 1551 ) can be solved analytically and matched 
to the exterior solution. The same procedure explained 
above then yields A/M 5 = — ||c -4 , which correctly re¬ 
produces the scaling, the sign, and also the order of mag¬ 
nitude of the Newtonian Love number in Eq. ( 1351 ) . This 
suggests that the different scaling of A and its opposite 
sign relative to the case of an ordinary fluid are due to 
the peculiar equation of state P(p) = —p of the de Sitter 
interior. 


We verified that the properties of the electric tidal Love 
numbers with £ > 2 are qualitatively similar to the £ = 2 
case discussed above, namely A ~ c~ 2e in the Newtonian 
limit and A ~ (1 — 2c) 2 in the BH limit for any £. 


Finally, the analogy with a Newtonian fluid having 
equation of state P(p) = — p is also useful to inter¬ 
pret the results obtained in the previous section for 
the quadrupole moment, namely that a slowly-spinning 
gravastar is prolate in the Newtonian limit whereas it 
becomes oblate approaching the BH limit]/]. For ordi¬ 
nary bodies, one can think of the oblateness as arising 
from the object’s response to the centrifugal pseudo¬ 
force. An object with negative pressure responds in 
the opposite direction that normal materials responds to 
forces (i.e., with an anti-restoring force directed as the 
external force). Since ordinary spinning objects become 
oblate, one might expect that an object with negative 
pressure will instead become prolate, in agreement with 
the Newtonian limit of Eq. m- Likewise, the fact that a 
gravastar becomes oblate at high compactness might be 
interpreted as due to the contribution of the thin shell 
(made of ordinary matter), which becomes dominant in 
this limit. 


We thank Leo Stein for suggesting this analogy. 
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FIG. 1. The I -A (top panel) and the Q-X (bottom panel) 
relations for a thin-shell gravastars as computed in this work 
and compared with those for a NS for various equations of 
state (see legend) spanning a wide range of NS deformability 
(cf. Refs. jlHi. [Ta| for details). The BH limit is denoted by a 
circular black marker. For NSs the curves are truncated at a 
value roughly corresponding to the maximum mass, whereas 
gravastar models connect continuously to the BH limit. Note 
that A < 0 in the gravastar case. 


IV. I-LOVE-Q RELATIONS FOR A GRAVASTAR 


Yagi and Yunes 00 discovered that for NSs and 
quark stars the dimensionless quantities 


I := 


I 

W’ 


Q 


Q 

X 2 M 3 ’ 



(36) 


satisfy nearly universal relations that are insensitive (at 
the percent level) to the equation of state of matter in 
the stellar interior. 

For the thin-shell gravastar model under consideration 
the dimensionless quantities (15S1) can be obtained analyt¬ 
ically as one-parameter functions depending only on the 
compactness [cf. Eqs. (0, (l2ll) and (1521) ]. These relations 
are valid in the full range 0 < c < 1/2, thus connecting 
smoothly the Newtonian regime to the BH limit. It is 
therefore interesting to compare our results with those 
obtained for NSs and BHs. Figure Q] shows a comparison 
of the I -A and the Q-X relations for thin-shell gravastars 
and for NSs with some realistic equations of state 00- 

We can observe various interesting features. First, 
while the relations for NSs are almost insensitive to the 
equation of state (on the scale of the plot the three curves 
lay on the top of each other), the same relations for 
gravastars are dramatically different in the entire range 
of compactness. We stress that our gravastar model cor¬ 
responds to a specific equation of state for the thin shell 
[cf. below Eq. (fljj) ] so that our results do not assess the 
universality (of absence thereof) of the I-Love-Q rela¬ 
tions for a thin-shell gravastar. Nonetheless, such results 


are sufficient to prove the opposite, i.e., that the I-Love- 
Q relations of a gravastar are different from those of a 
NS. This is due to the fact that for gravastars both the 
low-compactness and the high-compactness regimes are 
markedly different from the NS case. This is a nontrivial 
result, since for other compact stellar configurations (e.g. 
for quark stars) the I-Love-Q relations are the same as 
those for a NS for a variety of equations of state. 

Secondly, we stress that A is negative for gravastars, 
i.e. it has the opposite sign than for ordinary stars. As 
suggested by the Newtonian model discussed above, these 
peculiar properties are due to the exotic equation of state 
associated with a positive cosmological constant in the 
gravastar’s interior, P{p) = —p. 

Finally, in the gravastar case the mass quadrupole mo¬ 
ment does not have a definite sign, being negative for 
small compactness, becoming positive for c > 0.358 (or, 
equivalently, for | A| < 2.09), and approaching unity in the 
BH limit. Therefore, slowly-spinning gravastars are pro¬ 
late for moderately small compactness and become oblate 
as the compactness increases towards the BH limit. 
This behavior is qualitatively similar to that observed 
in strongly-anisotropic, incompressible NSs 0, 0. 


1. The approach to the BH limit 

As already discussed, the BH limit is continuous, 
namely 


/ —>■ 4, Q-> 1, Ah>0, (37) 


as c -> 1/2, precisely as in the Kerr case to quadratic 
order in the spin. Inverting Eq. ([6]), c = (2/) -1 / 3 , and 
using Eqs. (1^1) and (1M1) . we can easily compute the scal¬ 
ing of the I-X-Q relations near the BH limit, 


1 + 77 log(l-2(2/)- 1 / 3 
45 L 


A-s- — 


128 


( 2/)- 1/3 - £ 



(38) 

(39) 


as I —> 4 (or, equivalently, as c —> 1/2). In this specific 
model, the I-X-Q relations are nonpolynomial functions 
near the BH limit. This behavior is markedly different 
from what observed in strongly-anisotropic NSs. In that 
case, a given multipole moment A4 shows a critical be¬ 
havior near the BH limit, namely [0 0] 


<5 := ~Ta 1 ~ (1 — 2c) fc c ~ 1/2 , (40) 

AIbh 


where A4bh is the value of the multipole moment in the 
BH case and k is a critical exponent which depends only 
mildly (roughly within 10%) on the equation of state. For 
anisotropic NSs, the critical exponent for all multipole 
moments is roughly k ~ 4 0,0] . On the other hand, for 
gravastars we obtain k = 1 for the first current moment, 
Si := J = in, and we even obtain a nonpolynomial 

















behavior for the mass quadrupole moment, M 2 oc Q , for 
which 5 ~ [log(l — 2c)] -1 . 

Finally, we note that for gravastars /(c), Q(c ) and A(c) 
are monotonic functions of the compactness in the full 
range c £ [0,1/2], whereas they display a peculiar non¬ 
monotonic behavior in the case of strongly-anisotropic 

nss mm. 


2. Eccentricity profiles 


Recently, Ref. [32] proposed a phenomenological expla¬ 
nation for the nearly-universal relations of NSs in terms 
of an approximate symmetry emerging at high compact¬ 
ness. As the compactness increases, radial profiles of the 
eccentricity become nearly constant within the star, lead¬ 
ing to the emergence of isodensity self-similarity. 

It is therefore interesting to check whether the gravas- 
tar model enjoys a similar property. For our metric 
ansatz © and in the absence of a fluid, the eccentric¬ 
ity is simply defined as [Hi, IHoj] 

e(r) = V-3\M r ) -h(r)} ■ (41) 


For ordinary objects, V 2 — J 2 < 0 when the object is 
oblate, so that the eccentricity is a real number. In the 
case of a gravastar instead, V 2 — j -2 >0 in the interior 
the object for any compactness. Using the analytical re¬ 
sults previously presented, it is straightforward to obtain 
an analytical expression for the eccentricity in the entire 
space. The explicit form is not illuminating so we do 
not show it here, but in Fig. [5] we present the eccentric¬ 
ity profiles for various values of the compactness in the 
entire space. 

As anticipated, e 2 < 0 for r < R and any compact¬ 
ness. This is a peculiar feature, because it corresponds 
to an object which is (locally) prolate, even though the 
quadrupole moment can be either positive or negative 
depending on the compactness. Indeed, for c < 0.358 
(top panel of Fig. © we observe that e 2 < 0 in the en¬ 
tire space, in agreement with the fact that in this case 
Q < 0 and therefore the object appears prolate also at 
infinity. On the other hand, when c > 0.358 (bottom 
panel of Fig. © we observe that e 2 (r) changes sign in the 
exterior and becomes positive at large distances. In this 
case, the object has a positive quadrupole moment but a 
negative eccentricity in the interior. The turning point of 
the function e 2 (r) decreases as the compactness increases 
and tends to the radius as c —> 1/2. 

In the BH limit the eccentricity at the radius is dis¬ 
continuous. Indeed, for our gravastar model e 2 (R) —> 
— X 2 /4 as c —► 1/2, whereas the eccentricity of a (slowly- 
spinning) Kerr BH can be easily computed from Eq. (1411) 
and by using the external solution Col (CD with 5m = 
5q = 0. This yields 


2 3(r/M + l) 2 

e BH — 0 / /„r\4 X : 


( 42 ) 



r/R 


FIG. 2. Eccentricity profiles of a slowly-rotating, thin-shell 
gravastar. The top panel shows the profiles for prolate objects 
(Q < 0) which requires c < 0.358. In this case e 2 < 0 in the 
entire space. The bottom panel shows the profiles for oblate 
objects (Q > 0). In this case e 2 < 0 up to a certain point in 
the exterior, and then turns positive for larger values. The 
BH case, Eq. (1421) . is shown as a dashed orange curve. Note 
that the BH limit is discontinuous because e 2 (R) —> — x 2 /4 
as c —> 1/2, whereas e| H (r = 2 M) = 9x 2 /128. 


and the eccentricity at the event horizon is positive, 
e 2 (2 M) = 9x 2 /128. 

In light of these peculiar properties, it is difficult to 
make a connection to the eccentricity profiles of ordinary 
stars [26j, |32l , in which e 2 > 0 and the eccentricity re¬ 
mains nearly constant inside the star as the compactness 
increases. In the gravastar case, the eccentricity vanishes 
near the center, e 2 ~ r 2 as r —> 0, and does not approach 
a nearly constant profile as c —> 1/2, although it shows 
some interesting features that might be worth exploring 
in more details. In particular, an analysis of the eccen¬ 
tricity profiles for various equations of state of the thin 
shell is underway [53j . 


V. DISCUSSION 


The fact that the I-X-Q relations for gravastars are 
completely different from those of NSs and quark stars 
is not surprising. A gravastar is a rather exotic model 
whose interior is not described by an ordinary fluid. In 
particular, being a gravastar a purely gravitational ob¬ 
ject, rotation and tidal fields do not induce any fluid mo¬ 
tion. In fact, gravastars have been proposed to mimic 
BHs rather than compact stars and they share more fea¬ 
tures with the former, for example in the static case they 
can have any mass and a compactness c —> 1/2. 

More interestingly, the I-X-Q relations for gravas¬ 
tars approach smoothly the corresponding value in 
the (slowly-rotating) Kerr case as the compactness ap- 


8 (r/M) 4 
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proaches its maximum valucQ, c —> 1/2. This fact is re¬ 
markable because the BH limit does not need to be con¬ 
tinuous and, in fact, the eccentricity at the surface does 
not connect continuously to its Kerr value as c —> 1/2. 

The fact that the multipole moments approach their 
Kerr values has also some phenomenological consequence. 
It is usually assumed that a measurement of the mass 
quadrupolc moment of a supermassive compact object 
would be an ultimate test of the BH no-hair theorem 
and, in turn, of the “Kerr hypothesis” (cf. e.g. Ref. 
Gravastars challenge this paradigm, because they can be 
as massive and as compact as BHs and, furthermore, they 
have the same moment of inertia, mass quadrupole mo¬ 
ment and even the same (vanishing) tidal Love numbers. 

Although we did not compute higher multipole mo¬ 
ments and the magnetic Love numbers, there is no reason 
to expect qualitatively different results in those cases. It 
is therefore natural to conjecture that the exterior geom¬ 
etry of a spinning, thin-shell gravastar is equivalent to 
that of a Kerr BH in the c —> 1/2 limit. This would 
imply that gravastars and BHs are indistinguishable on 
the basis of tests probing the multipolar structure of the 
geometry (e.g. through geodesic motion) or probing the 
tidal deformability of these objects. Our results suggest 
that the only direct waj0 to distinguish these objects 
from BHs is through gravitational-wave observations of 
the ringdown [40l 4H, 43 or of the inspiral in the extreme 
mass ratio limit [59l . Elil- 

On the other hand - if one is willing to consider the 
gravastar proposal seriously - the same formation mech¬ 
anism of BH-like gravastars should also produce stellar- 
like objects with smaller compactness and with a mass 
comparable to that of a NS. The existence of different I- 
X-Q relations can therefore be used to rule out stellar-like 
gravastars once independent measurements of any two 
elements of the triad become available. Figure [T] shows 
that, for a fixed value of |A|, the deviations of I and Q for 
a gravastar relative to a NS are larger than ~ 200% and 
~ 80%, respectively. For moderately compact gravastars, 
the dimensionless moment of inertia can be larger than 
that of a NS with same quadrupole moment by more 
than a factor 10. Therefore, even a measurement with 
large uncertainties would easily discriminate between a 


gravastar and a NS. 

Finally, our results provide a simple testbed to ex¬ 
plore the approach of the I-X-Q relations to the BH 
limit. In this limit the behavior of gravastars is remark¬ 
ably different from that of anisotropic NSs discussed in 
Refs. [25], (26j. In particular, in the BH limit gravastars 
do not satisfy the critical behavior found for anisotropic 
stars, either by displaying a different critical exponent 
or even by displaying a nonpolynomial behavior near the 
critical point. Furthermore, the eccentricity inside the 
star does not approach a nearly-constant profile as the 
compactness increases, but displays some peculiar fea¬ 
tures which are absent in the NS case [26l. |32|. 

It would be interesting to extend this analysis to other 
“BH mimickers”, e.g. wormholes, quasi-BHs, other black 
foils (cf. e.g. Refs. [fHHU) and to gravastars with dif¬ 
ferent equations of state for the thin shell, to understand 
whether the approach to the BH limit is characterized by 
some universal behavior. 

Note added: After a first version of this work was sub¬ 
mitted, we became aware of a very similar analysis in 
Ref. |52[, which also investigates slowly-rotating thin- 
shell gravastars to quadratic order in the spin. Our work 
and that of Ref. [52| are complementary to each other: 
we investigate the I-X-Q relations and the approach to 
the BH limit, whereas Ref. [ 52 } is devoted to a rigorous 
analysis of the junction conditions at the interface and 
on the properties of the thin shell. 
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